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In  this  paper  a  dynamic  model  of  a  single  solid  oxide  fuel  cell  (SOFC)  is  developed  using  a  volume  element 
methodology.  It  consists  of  a  set  of  algebraic  and  ordinary  differential  equations  derived  from  physical 
laws  (e.g.,  the  first  law  of  thermodynamics,  Fields  law,  and  Fourier’s  law),  which  allow  for  the  prediction 
of  the  temperature  and  pressure  spatial  distribution  inside  the  single  SOFC,  as  functions  of  geometric 
and  operating  parameters.  The  thermodynamic  model  is  coupled  with  an  electrochemical  model  that 
is  capable  of  determining  the  voltage,  current,  and  power  output.  Based  on  the  simulation  results,  the 
internal  configuration  (structure  of  the  positive  electrode-electrolyte-negative  electrode  assembly)  and 
the  operating  conditions  (air  stoichiometric  ratio  and  fuel  utilization  factor),  as  well  as  their  impact  on  the 
performance  of  the  single  SOFC  are  discussed.  Optimal  geometric  and  operating  parameters  are  obtained 
so  that  electrical  power  of  the  single  SOFC  at  the  nominal  operating  point  is  maximized.  The  method 
used  is  general  and  the  fundamental  optimization  results  are  sharp,  showing  up  to  a  357%  single  SOFC 
performance  variation  within  the  studied  parameters’  range,  therefore  these  findings  show  the  potential 
to  use  the  model  as  a  tool  for  future  SOFC  design,  simulation  and  optimization. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  electrochemical  devices  that  convert  chemical 
energy  in  fuels  directly  to  electricity.  Unlike  the  conventional  power 
generation  methods,  there  are  no  intermediate  steps  such  as  pro¬ 
ducing  heat  by  burning  fuels  and  converting  it  to  mechanical 
work  with  a  heat  engine.  Thus,  the  fuel  cells  are  not  restricted 
by  some  thermodynamic  limitations  and  have  high  energy  effi¬ 
ciency.  Besides,  except  for  the  small  amount  of  C02  emitted  as 
the  byproduct  of  hydrogen  generation  from  hydrocarbon  com¬ 
pounds,  fuel  cells  exert  minimum  impact  on  the  environment. 
Because  of  all  these  virtues,  they  are  considered  as  one  of  the 
most  attractive  solutions  for  clean  power  systems  and  receive 
widespread  attention  from  the  automotive  and  power  genera¬ 
tion  industry,  as  well  as  many  research/academic  organizations 
[1]. 

Among  different  types  of  fuel  cells,  solid  oxide  fuel  cells  (SOFCs) 
are  particularly  attractive  due  to  their  special  features:  great  fuel 
flexibility,  the  feasibility  of  fuel  internal  reforming,  high  energy 
conversion  efficiency,  low  emission  and  noise,  and  the  potential 
application  in  the  highly  efficient  cogeneration  power  systems. 
However,  there  are  key  challenges  to  be  overcome  before  their  mass 
production  and  widespread  application,  such  as  power  density  aug¬ 
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mentation,  operating  condition  control  and  fuel  cell  configuration 
optimization. 

Much  research  and  development  has  been  devoted  to  SOFCs  and 
published  in  the  scientific  literature.  Virkar  et  al.  [2]  discussed  the 
effect  of  electrode  microstructure  on  activation  and  concentration 
potential  losses  in  the  SOFC,  where  the  concept  of  effective  charge 
transfer  resistance  was  introduced.  Chan  et  al.  [3]  presented  a  com¬ 
plete  polarization  model  of  SOFC  that  could  work  under  different 
design  and  operating  conditions.  Sorrentino  et  al.  [4]  built  a  ID 
model  for  the  planar  SOFC  and  validated  it  with  the  experimen¬ 
tal  data  from  Pacific  Northwest  National  Laboratory.  In  the  paper, 
the  effects  of  the  several  important  operating  variables  on  SOFC 
performance  were  evaluated,  such  as  the  fuel  utilization  factor, 
operating  temperature  and  pressure.  Vargas  et  al.  [5]  built  a  ID 
model  for  a  PEM  fuel  cell,  and  with  a  structured  procedure,  opti¬ 
mized  the  internal  and  external  structure  of  the  PEM  fuel  cell  to 
achieve  the  maximum  net  power  under  different  stoichiometric 
ratios.  The  model  was  further  validated  with  the  experiment  data 
by  Martins  et  al.  [6].  Based  on  the  model  of  the  unit  PEM  fuel  cell, 
Vargas  et  al.  [7]  developed  a  model  for  PEM  fuel  cell  stacks  and  the 
effect  of  internal  and  external  configuration  on  stack  performance 
was  discussed.  For  SOFCs,  Ordonez  et  al.  [8]  optimized  geometric 
features  of  a  single  SOFC  based  on  a  simplified  model. 

Young  [9]  and  Ma  et  al.  [10]  recently  reviewed  the  status  of 
three-dimensional  SOFC  modeling.  Young  [9]  discussed  current 
thermo-fluid  modeling  capabilities  for  PEM  fuel  cell  and  SOFCs, 
and  highlighted  key  challenges  such  as  gas  diffusion  in  the  porous 
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Nomenclature 

A  area,  m2 

A' 0pt  area  of  the  active  TPB  regions,  m2 

c  molar  specific  heat,  J  mol-1  K-1 

cv  specific  heat  capacity  at  constant  volume, 

J  mol-1 1<-1 

cp  specific  heat  capacity  at  constant  pressure, 

Jkg-’K-1 

D  mass  diffusivity,  m2  s-1 

Dh  gas  channel  hydraulic  diameter,  m 

Dp  average  diameter  of  the  pore,  m 

Dw  width  of  a  corrugation,  m 

E  activation  energy,  J  mol-1 

/  friction  factor 

F  Faraday’s  constant,  Cequiv.-1 

h  heat  transfer  coefficient,  W  m-2  K-1 

H  enthalpy,  J  mol-1 

i  current  density,  A  nrr 2 

i'o  exchange  current  density,  A  m-2 

/as  limiting  current  density(anode-limited),  A  m-2 

ics  limiting  current  density(cathode-limited),  Am-2 

I  working  current,  A 

k  thermal  conductivity,  W  m- 1 1<-1 

L  length,  m 

U  reaction  region  thickness,  m 

L' opt  active  TPB  region  thickness,  m 

m  mass,  kg 

fn  mass  flow  rate,  kg  s-1 

M  molecular  weight,  kg  mol-1 

n  equivalent  electron  per  mole  of  reactant, 

equiv.  mol-1 

fl  molar  rate  involved  in  the  electrochemical  reaction, 

mols-1 

flin  molar  flow  rate  at  the  inlet,  mol  s-1 

n0ut  molar  flow  rate  at  the  outlet,  mol  s-1 

nCh  number  of  channels 

nm  i  number  of  moles  of  component  i  in  the  mixture 
nm, t  total  number  of  moles  in  the  mixture 

N  heat  capacity,  J  I<-1 

p  pressure,  Pa 

P  dimensionless  pressure 

Pr  Prandtl  number 

Q  heat  transfer  rate,  W 

R  universal  gas  constant,  J  mol-1 1<-1 

Rj  specific  gas  constant,  J  kg-1  K-1 

Re  Reynolds  number 

S°  standard  entropy,  J  mol-1 1<-1 

t  time,  s 

T  temperature,  K 

u  velocity,  ms-1 

V  voltage,  V 

Vj  total  volume  of  the  single  SOFC,  m3 

W  power,  W 

We  electric  power,  W 

WN  net  power,  W 

Wp  pumping  power,  W 

X  molar  fraction  of  component  in  the  mixture 

Zr  grain  size,  m 

Greek  symbols 

a  fraction  of  the  reaction  heat  that  is  generated  at  the 

anode 

ac  charge  transfer  coefficient 


at  anode  thermal  diffusivity 

P  ohmic  resistance,  £2 

X  gas  channel  aspect  ratio 

8  pre-exponential  factor,  nrr2 

s  emissivity 

0  porosity 

Yf  fuel  utilization  factor 

p  overpotential,  V 

0  dimensionless  temperature 

§  dimensionless  length 

r  dimensionless  time 

Xair  stoichiometric  ratio 

/x  viscosity,  Pa  s 

v  diffusion  volume  of  simple  molecules,  cm3 

p  density,  kg  nr 3 

a  electrical  conductivity,  m-1 

orrad  Stefan-Boltzmann  constant 

0  tortuosity 

k  isentropic  index 

Subscript  and  superscript 
0  initial  condition 

1  bipolar  plate  (interconnect) 

2  fuel  channel  or  fuel 

3  anode 

4  electrolyte 

5  cathode 

6  air  channel  or  air 

~  dimensionless 

a  anode 

act  activation 

actconc  activation  and  concentration 

active  the  active  reaction  region 

air  air 

c  cathode 

ch  channel 

end  conduction 

env  convection 

cone  concentration 

eff  effective 

f  fuel 

h  hydraulic 

in  inlet 

H2  hydrogen 

H20  water 

in  inlet 

limit  limit 

mm  two-way  maximum 

mmm  three-way  maximum 

mass  mass  transfer 

N2  nitrogen 

n  nominal  operating  point 

o  standard  condition  [gases  at  1  atm,  298.1 5  K] 

02  oxygen 

ocp  open  circuit  potential 

ohm  ohmic 

opt  optimal 

out  outlet 

rad  radiation 

ref  reference 

ret  reaction 

t  channel  shoulder 

total  total 
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TPB  three-phase  boundary 

x  x  direction 

y  y  direction 

z  z  direction 


layer,  the  internal  fuel  reforming  rate  and  the  transport  of  gas 
species,  ions,  and  electrons  accompanied  with  the  electrochem¬ 
ical  reaction.  Ma  et  al.  [10]  reviewed  the  current  computational 
fluid  dynamics  (CFD)  modeling  of  fuel  cells.  Bove  and  Ubertini  [11] 
developed  a  complete  three-dimensional,  time-dependent  SOFC 
numerical  model  which  accounts  for  various  phenomena  occurring 
in  each  component  of  the  fuel  cell.  Janardhanan  and  Deutschmann 
[12]  made  CFD  analysis  of  an  internally  reforming  anode  sup¬ 
ported  SOFC  button  cell,  which  accounts  for  the  heterogeneous  fuel 
reforming  chemistry,  electrochemistry  and  porous  media  trans¬ 
port  inside  the  fuel  cell.  In  the  paper  the  effect  of  steam  content 
in  the  anode  feed  stream  on  the  resulting  overpotential  losses  and 
surface  coverages  of  various  species  at  the  three-phase  boundary 
were  discussed.  Furthermore,  Janardhanan  et  al.  [13]  introduced  a 
mathematical  model  to  calculate  the  volume  specific  three-phase 
boundary  length  in  the  porous  composite  electrodes  of  SOFC,  which 
is  of  crucial  importance  to  the  SOFC  performance. 

Though  very  useful  refined  spatial  information  can  be  extracted 
from  three-dimensional  models,  the  computational  time  require¬ 
ments  prevent  their  use  for  optimization  strategies  in  which 
multiple  runs  for  different  geometries  and  operating  conditions  are 
necessary.  As  a  result,  there  is  a  need  for  an  efficient  computational 
methodology,  which  captures  the  dominant  physical  and  electro¬ 
chemical  phenomena  occurring  in  SOFC  and  is  accurate  enough 
to  depict  the  response  of  the  single  SOFC  under  different  internal 
configurations  and  operating  conditions. 

In  this  paper  a  one-dimensional  single  SOFC  model  based  on  the 
fundamental  physical  laws  is  presented.  The  single  SOFC  is  divided 
into  six  lumped  control  volumes  that  correspond  to  the  most  repre¬ 
sentative  parts  of  a  single  fuel  cell.  The  phenomena  such  as  heat  and 
mass  transfer,  fluid  flow,  electrochemical  reaction,  potential  losses 
and  power  consumption  and  generation  are  taken  into  account, 
and  the  solution  includes  the  temperature  and  pressure  field  in 
each  control  volume,  the  polarization  curves  and  the  power  output 
of  the  single  SOFC.  Based  on  the  simulation  results  obtained  with 
the  model,  the  effects  of  operating  parameters,  such  as  air  stoichio¬ 
metric  ratio  and  fuel  utilization  factor,  as  well  as  the  PEN  (positive 
electrode-electrolyte-negative  electrode  assembly)  internal  struc¬ 
ture  on  SOFC  performance  are  investigated.  Finally,  the  optimal 
values  for  each  investigated  design  and  operating  parameter  of  the 
single  SOFC  are  proposed  for  maximum  fuel  cell  performance. 

2.  Thermal  model 

The  configuration  of  a  single  SOFC  is  shown  in  Fig.  1.  It  consists 
of  six  control  volumes:  the  solid  bipolar  plate— interconnect  (CV1 ), 
the  fuel  channels  (CV2),  the  anode  layer  (CV3),  the  electrolyte  layer 
(CV4),  the  cathode  layer  (CV5),  and  the  oxidant  channels  (CV6). 

The  fuel  and  oxidant  supplied  at  the  channels  diffuse  through 
the  anode  layer  and  cathode  layer,  respectively,  where  they  react 
separately  at  the  ionic  conductor/electronic  conductor/gas  inter¬ 
face  (three  phase  boundaries,  TPB).  The  bipolar  plate  (interconnect) 
provides  a  series  of  electrical  connectors  to  the  adjacent  cells  or 
to  the  external  circuit,  and  at  the  same  time  serves  as  a  gas  bar¬ 
rier  between  the  fuel  channels  and  oxidant  channels  of  adjacent 
cells. 

A  wide  range  of  fuels  could  be  used  in  the  electrochemical  reac¬ 
tion  of  SOFC,  including  various  hydrocarbon  fuels.  For  simplicity,  in 


fuel 


Fig.  1.  The  3D  configuration  of  a  planar  single  SOFC. 


this  paper  pure  hydrogen  and  air  are  considered  as  fuel  and  oxidant, 
respectively. 

Due  to  the  heat  and  mass  interaction  inside  the  single  SOFC,  the 
thermal  model  includes  mass  and  energy  conservation  equations 
for  each  control  volume,  as  well  as  the  electrochemical  reaction 
equations  at  the  two  electrodes.  The  electrochemical  model  uses 
the  temperature  and  pressure  fields  in  the  evaluation  of  the  per¬ 
formance  of  the  single  SOFC. 

The  electrochemical  model  contains  the  analysis  and  computa¬ 
tion  of  reversible  electrical  potential  and  various  kinds  of  losses, 
such  as  the  losses  due  to  surface  overpotentials  (poor  electrocatal¬ 
ysis),  slow  diffusion  and  all  internal  ohmic  losses  through  the  cell. 
The  actual  cell  potential  results  from  the  reversible  electrical  poten¬ 
tial  by  subtracting  all  the  losses,  which  are  functions  of  the  working 
current,  temperature  and  pressure,  and  SOFC  configuration.  The 
computation  of  the  actual  electrical  power  output  from  the  SOFC 
and  the  power  consumed  pumping  fuel  and  oxidant  into  the  gas 
channels  are  also  included  in  the  model  to  obtain  the  net  power 
output. 

2A.  Mass  balances  in  the  fuel  and  air  channels 

In  a  single  SOFC,  there  are  electrochemical  reactions  both  in  the 
anode  layer  and  the  cathode  layer.  These  reactions  are  expressed 
by  the  equations: 


H20  +  2e  (anode) 

(i) 

20=  (cathode) 

(2) 

The  hydrogen  in  the  anode  layer  and  the  oxygen  in  the  cath¬ 
ode  layer  are  supplied  by  the  fuel  channels  and  oxidant  channels, 
respectively.  Based  on  the  electrochemical  reaction  equations,  the 
molar  flow  rate  of  hydrogen  at  the  inlet  and  outlet  of  fuel  channels, 
as  well  as  hydrogen  consumed  by  the  electrochemical  reaction,  are 
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determined  by  the  working  current  (/): 

II  I  (I  A 

n^’in  -  -  Wf’  nH2’out  ~  Wf  ~  ) 

I  I 


I 

2 F 


0) 


where  n  is  the  number  of  mols  of  electrons  transferred  per  mol  of 
reactant,  for  the  hydrogen  n  =  2,  and  for  the  oxygen  n  =  4.  F  is  the 
Faraday  constant,  96500  C  mol-1 ,  and  yf  (yf  <  1 )  is  the  fuel  utiliza¬ 
tion  factor  and  defined  as  the  ratio  between  the  fuel  consumed  in 
the  reaction  and  the  fuel  supplied  in  the  fuel  channels. 

In  the  fuel  channels,  the  average  molar  fractions  Xi  =  nm<ilnmj, 
defined  as  the  average  number  of  mols  of  component  i  divided  by 
the  total  number  of  mols  of  the  mixture,  is  approximated  as: 


y  Yf  y  -i  Yf 

ah2o  =  ~2’  Ah2  =  1  ~  ~2 


(4) 


Similarly,  the  molar  flow  rate  of  oxygen  at  the  inlet  and  outlet 
of  air  channels  and  the  oxygen  consumed  by  the  electrochemical 
reaction  are: 


/  .  /  . 

n02,in  -  npAair  -  air > 


I 


n02  =  = 


I 


nF  4  F 


^02,out  —  (^air  1)^-  (^air  1)^ p  ’ 

(5) 


where  Aair  is  the  air  stoichiometric  ratio,  and  defined  as  the  air 
provided  in  the  channels  divided  by  the  air  needed  for  the  electro¬ 
chemical  reaction  of  interest. 

In  the  present  study,  the  oxidant  provided  in  the  oxidant  chan¬ 
nels  is  air,  which  is  assumed  to  be  composed  of  oxygen  and  nitrogen 
only.  The  molar  flow  rate  of  nitrogen  at  the  inlet  and  outlet  are 
obtained  as: 


/Aair^N2,in  ^air^N2,in 
J2,in  =  nN2.out  =  nFv  -  =  4Fv  - 


(6) 


The  average  molar  fractions  of  oxygen  and  nitrogen  in  the  oxi¬ 
dant  channels  are 


*02  = 


2Aair^o2,in  ^02,in  ^o2,in 


*n2  = 


2(^air  ^02,in) 

2A-air^N2,in  —  ^02,in^N2,in 
^(•^air  —  ^02,in) 


(7) 


2.2.  Energy  balances 
2.2.1.  Dimensionless  variables 

The  application  of  the  first  law  of  thermodynamics  to  every  con¬ 
trol  volume  of  the  single  SOFC  illustrated  in  Fig.  1  yields  expressions 
for  their  temperature  variations  during  the  heat  and  mass  trans¬ 
fer  processes  associated  with  the  SOFC  operation.  In  addition,  to 
simplify  the  acquired  equations,  they  are  written  in  dimensionless 
form  and  all  the  quantities  in  them  are  rescaled  to  values  as  close  as 
possible  to  1,  which  in  turn  simplifies  the  numerical  calculations. 

In  a  dimensionless  model,  all  variables  are  directly  proportional 
to  the  actual  dimensional  ones.  Therefore,  this  allows  for  scaling 
up  or  down  any  system  with  similar  characteristics  to  the  system 
analyzed  by  the  model.  Another  important  aspect  is  that  any  dimen¬ 
sionless  variable  value  used  in  the  simulations  could  represent  an 
entire  and  numerous  set  of  dimensional  values  by  varying  appro¬ 
priately  the  parameters  in  the  dimensionless  variables  definition, 
which  by  itself  stresses  the  generality  of  the  dimensionless  model. 
Physically,  the  set  of  results  of  a  dimensionless  model  actually 
represents  the  expected  system  response  to  numerous  combina¬ 
tions  of  system  parameters  (geometry,  architecture)  and  operating 


conditions  (e.g.,  ambient  conditions  and  mass  flow  rates),  with¬ 
out  having  to  simulate  each  of  them  individually,  as  a  dimensional 
model  would  require.  Based  on  that  reasoning,  it  is  convenient  to 
search  for  an  alternative  formulation  that  eliminates  the  physical 
dimensions  of  the  problem. 

For  instance,  all  the  actual  lengths  in  the  model  are  nondimen- 
sionalized  in  the  following  way,  from  which  their  corresponding 
dimensionless  lengths  are  obtained: 

,  i  =  1, 2,  3,  4,  5,  6,  t,  ch,  x,y,  z 

V. T1/J 

where  VT  =  LxLyLz  represents  the  total  volume  of  the  single  SOFC. 
Similarly,  the  dimensionless  temperatures  are  defined  as: 

0,  =  fi-,  i=  1,2,3, 4,  5,  6 

'ref 

where  Tref  is  referenced  to  the  ambient  temperature.  The  dimen¬ 
sionless  pressures  are: 

Pi=—,  i  =  H2,  02,  H20,  N2,  air,  f,  TPB 
P  ref 

where  pref  is  referenced  to  the  ambient  pressure. 

Other  dimensionless  quantities  include: 

The  dimensionless  time  r  =  t/tref,  where  tref  =  (4V^3)lat  and 
at  =  /<3/(P3  cp>3 )  refers  to  the  thermal  diffusivity  of  the  anode. 

The  dimensionless  areaAz-  =  Aj/(1/^3). 

The  dimensionless  mass  fri/  =  where  mref  =  p3VT. 

The  dimensionless  mass  flow  rate  rhj  =  m//mref,  where  mref  = 

(p3VT)/tref- 

The  dimensionless  thermal  conductivity  /q  =  (/qV^3)/(rhrefCP)f). 
The  dimensionless  convective  heat  transfer  coefficient  = 

(h,\/T2/3)/(mrefcp,f). 

The  dimensionless  Stefan-Boltzmann  constant  arad  = 

(^radVT2/3rr3ef)/(mrefCp,f). 

The  dimensionless  heat  transfer  rate  and  power  Q  = 
Qil(mrefCp  fTref)  and  Wj  =  14^7(rhrefCp  fTref). 

In  all  the  equations  above,  the  subscript  i  indicates  a  substance 
or  a  location  in  the  fuel  cell. 


2.2.2.  Energy  conservation 

The  thermal  model  consists  of  six  ordinary  differential  equations 
which  are  derived  by  applying  the  energy  conservation  law  to  each 
control  volume.  In  the  present  study  the  single  SOFC  is  considered 
to  be  an  internal  cell  in  a  SOFC  stack  and  has  the  interaction  with  its 
adjacent  cells.  It  is  also  assumed  that  the  adjacent  cells  have  similar 
internal  temperature  profiles. 

For  CV1,  energy  conservation  states  that: 


d0\  _  (Qrad,l— 3  +  Qrad,l-5  +  Qcnv,l-2  +  Qcnv,l-6  +  Qohm,l  +  Qcnd,l-3  +  Qcnd,l-5) 


dr 


Ni 


(8) 


where  the  dimensionless  heat  capacities,  Nb  are  given  as: 


Ni  = 


wiiCpi 

mref^p,f 


i=  1,3,  4,  5 


(9) 


where  m\  =  p\  [L\LyLz  +  ( lyl2  +  L6Lz)Lt(nCh  + 1 )]  and  m*  =  PiLyLzLit 
1  =  3,4,  5. 

In  Eq.  (8),  the  radiation  heat  transfer  from  the  anode  layer  and 
cathode  layer  to  the  bipolar  plate  are: 


Qrad.W  =  ^rad AadW0/  “  J  =  3>  5  (10) 

where  ,Arad  3  =  (nchLchIy)/(V2/3)  and/tradj5  =  (nchIchLz)/(V,2/3)  are 
the  dimensionless  radiation  heat  transfer  areas  between  the  elec¬ 
trodes  and  the  bipolar  plate. 
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The  effective  emissivities  are  estimated  by  [14]: 

1 

6eff.3  -  i/fi3  +  (Lch/(Lch  +  2L2))(l/fi,  -  1) 

(11a) 

1 

Seff’5  “  l/e5  +  (WUci,  +  2L6))(1/s1  -  1) 

(lib) 

In  Eq.  (8),  the  heat  conduction  from  the  anode  layer  and  the 
cathode  layer  to  the  bipolar  plate  are: 

a  2Acnd  3  (03  -  0\  )/<3/<i 

(2cnd,l-3  =  ,,  r  ,  r  . 

(§3/<l  +  §2/<3) 

(12a) 

a  2Acnd  5(05  -  0\  )/<5/<  1 

(2cnd,l-5  =  ,,  r  ,  r  . 

(§5^1  +  ^<5) 

(12b) 

where  Acnd  3  —  (1  03)(^ch  +  1  )£t Hy  ^tld  ^cnd,5  —  0  05)(^ch  T 
1  )§t Hz  stand  for  the  conduction  heat  transfer  areas  between  bipolar 
plate  and  electrodes.  </>3  and  </>5  represent  the  anode  and  cathode 
porosities,  respectively. 

In  addition,  the  dimensionless  convective  heat  transfer  from  fuel 
channels  and  oxidant  channels  to  the  bipolar  plate  are  expressed 
as: 

Qcnv, i  —j  =  hjAcnv  j{0j  -  ),  j  =  2,6  (13) 

where  Acnv,2  =  nch(£ch  +  2§2 )£y  and  Acnv,6  =  nch(?ch  +  2§6 )&• 

The  convective  heat  transfer  coefficients  hj  in  the  fuel  and  oxi¬ 
dant  channels  are  estimated  according  to  the  flow  regime.  For  the 
laminar  flow  {Reh  <  2300)  [  1 5 ] : 

fjRehJ  =  24(1  -  1.3553Xj  +  1.9467x,2  -  1.7012x3  +  0.9564x4 
-0.2537xf),  j  =  2,6  (14a) 


In  the  present  study,  the  fuel  channels  are  open  control  vol¬ 
umes  filled  with  hydrogen  and  its  oxidation  product,  i.e.,  water.  It 
is  also  assumed  that  the  total  pressure  in  the  fuel  channels  and 
oxidant  channels  are  known  and  constant.  Applying  the  energy 
conservation  law  to  CV2  yields: 

d02  _  [— Qcnv, 1-2  +  Qcnv, 2-3  +  0mass,2 ]  ('ll) 

~dr  ~  W2 


where  the  dimensionless  heat  capacity  of  CV2,  N2,  is  given  as  fol¬ 
lows: 


~  _  PfnchLchL2Ly(XH2Cv,H2  +^h2oCv,h2o) 
RT2mrefCp  f 


(18) 


Qcnv, 2-3  is  the  dimensionless  convective  heat  transfer  between 
the  anode  and  the  fuel  channels  and  calculated  by:  Qcnv, 2-3  = 
h2nCh§ch?y(^  -  03)(^3  -  ^2).  where  the  convective  heat  transfer 
coefficient,  h2,  is  computed  with  Eq.  (14).  Qmass,2  stands  for  the 
energy  variation  associated  with  the  mass  transfer  (fuel  and  water 
vapor)  in  CV2  when  they  cross  the  boundaries: 


0mass,2 


”H2,in(^H2,f  -  ^H2 , 2 )  +  ^H2o(^H20,3  ~  ^H20,2) 
^ref^p,f^ref 


(19) 


where  hH2o  =  Hh2- 

Heat  is  generated  when  the  electrochemical  reactions  expressed 
by  Eqs.  (1 )  and  (2)  take  place.  In  the  present  study,  it  is  assumed  that 
the  total  thermal  effect  of  the  two  semi-reactions  on  the  tempera¬ 
ture  profile  of  the  single  SOFC  can  be  approximated  as  one  overall 
electrochemical  reaction:  H2  +  (l/2)02  H20,  which  takes  place  at 

the  electrolyte  temperature  T4.  The  total  heat  generation  is  given 
by  [17]: 


=  7.541(1  -  2.610Xj  +4.970Xj2  -  5.1 19xf  +  2.702x4 


-0.548X3),  j  =  2,6 

(14b) 

For  the  turbulent  flow  (2300  <  Reh  <  2  x  1 04 )  [  1 6] : 

fj )  =  0.079  Relf,  j  =  2,6 

(14c) 

hjDhj  (fj/2)(Rehj  -  1000)PTj  . 

ki  1  +  12.7(/j/2)1/2(Pr2/3  -  1)'  1 

(14d) 

where  X/  =  FCh/Fj,  when  Lch<Lj  and  Xj  =  Lj/Lc h*  when  Lch>Lj; 
Dhj  =  2LchLj/(Lch  +  Lj)\  Pr  is  the  gas  Prandtl  number;  /is  the  friction 
factor. 

The  Reynolds  number  is  a  function  of  the  average  flow  velocities 
in  the  channels  Uj ,  which  are  approximated  in  terms  of  the  molar 
flow  rates: 

fi[^H2,in7f,in  +  (^H2,out  +  ^H20,out)72  ] 

2PfTChT2nch 

(15a) 

R[(flo2,in  +  ^N2,in)Zair,in  +  (^02,out  +  ^N2,out)^6] 
2Pair^ch^6^ch 

(15b) 

The  dimensionless  ohmic,  activation  and  concentration  heating 
rates  are  given  as  follows: 

QohmJ=^°,hmi  l  j=h  3,4,5 

mref^p,f^ref 

(16a) 

n  ^act’J  +  Pconc,j  T  .  Q  £- 

CiactconcJ  —  ■  T  J  ~  D 

^ref^p,flref 

(16b) 

where  the  activation,  concentration  and  ohmic  potential  losses  in 
each  compartment,  (pactj,  pconcj-,  p0 hmj),  are  discussed  later  in  the 
electrochemical  model. 


Orel  = 


^H2[-^4  (Sh2qi4  _SH2,4  “  (V2)SS2,4)  +RT4  ln(PH20,TPB/PH2,TPB(Po2JPB)  7  )] 


ft1refCp,fT-ef 


(20) 


Part  of  the  dimensionless  reaction  heat  rate,  aOrct,  is  generated 
at  the  anode  (0  <  a  <  1 ),  and  the  rest,  ( 1  -  oOQrct,  is  generated  at  the 
cathode. 

According  to  the  assumption,  the  standard  entropy  of  forma¬ 
tion  S°  in  Eq.  (20),  is  evaluated  at  the  electrolyte  temperature 
T4.  The  three-phase  boundaries  pressures,  Ptpb,  refer  to  the  par¬ 
tial  pressures  of  reactants  and  products  at  the  place  where  the 
electrochemical  reaction  actually  takes  place.  Equations  for  the 
three-phase  boundaries  pressures  are  given  by  Kim  et  al.  [18]  and 
Aguiar  et  al.  [  1 9]  with  the  assumption  that  three-phase  boundaries 
pressures  are  uniform  within  the  reaction  regions  and  approxi¬ 
mated  as  the  partial  pressure  at  the  electrode/electrolyte  interface: 

n  n  Y  I  RT3^3  LsI 

Ph2o,tpb  =  Pf*H2o  +  2 FD^y7z  (21a) 


Ph2,TPB  =  Pf^H2  - 


RT 3^3  L3/ 
2FD303  LyLz 


(21b) 


Po2,TPB  —  Pair  “  (Pair 


^02Pair)exp 


(  RT5xlf5L5I  \ 

\4FD505pairLyLz  J 


(21c) 


where  ^  is  the  tortuosity,  D  is  the  diffusivity,  pf  is  the  total  pressure 
in  the  fuel  channels,  pair  is  the  total  pressure  in  the  oxidant  channels. 
According  to  Kim  et  al.  [18],  binary  diffusion  dominates  the  gas 
transport  through  the  porous  electrodes.  The  diffusivities,  Dj,  in  the 
anode  layer  and  the  cathode  layer  are  obtained  from  [20]: 


1.43  x  10~7Tj  75(Mh2  +Mh2q)1/2 

P f(  2JWh2  Mh2q  )’ /2  ( vfj  +  <30)2 


(22a) 
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Fig.  2.  Model  of  the  positive  electrode-electrolyte-negative  electrode  assembly  (PEN)  internal  structure  with  the  composite  electrodes. 


1.43  X  +Mn2)1/2 

Pair(2M02MN2)1/2(vf23  +  <3)2 


(22b) 


where  v  is  the  diffusion  volume  of  simple  molecules. 

Therefore,  for  CV3,  the  energy  balance  accounts  for  energy 
change  due  to  the  heat  transfer,  potential  losses  heating,  mass 
transfer,  and  the  heat  generated  during  the  electrochemical  reac¬ 
tion,  as  follows: 


The  energy  change  due  to  mass  transfer  across  the  electrolyte 
boundaries  is  calculated  as: 


Qm, 


iass,4 


fIH2(^H2,3  -  ^H2,4)  +  iio2(^02,5  -  ^02,4) 
^ref£p,f^ref 


(28) 


Similarly  to  the  anode  layer,  the  transient  variation  of  the 
cathode  layer  temperature  is  obtained  based  on  the  first  law  of 
thermodynamics  as  follows: 


dOs 

dr 


[— Qrad,l-3  —  Qcnv,2-3  +  Qohm,3  +  Qactconc,3  +  Qcnd,3-4  _  Qcnd,l-3  +  Qmass,3  +  ^Orct] 


n3 


(23) 


d0s 

dr 


[_ Qrad,l-5+Qcnv,5-6+Qohm,5+Qactconc,5— Qcnd,4-5_ Qcnd,l-5+Qmass,5  +  0  —  C*0Qrct] 

_ 


(29) 


where  Qcnd,3-4  is  the  dimensionless  heat  conduction  from  the  elec¬ 
trolyte  to  the  anode  layer,  and  it  is  given  by: 


Qcnd,3-4  — 


2k3MCnd,3— 4(04  ~  03 ) 

£4/(3  +  ^3^4 


(24) 


The  effective  heat  conduction  area,  Acnd  3_4,  in  Eq.  (24)  is  given  by. 
A:nd,3-4  =  0  -  03  )§y?z  Besides,  Qmass,3,  the  energy  change  of  CV3 
due  to  the  mass  flow  across  its  boundaries,  is  given  by: 


Qmass,3  — 


hH2(Hn2,2  -  ^H2,3 )  +  hH2o(^H20,4  -  #H20,3 
hlref^p,f^ref 


(25) 


The  energy  equation  for  CV4,  based  on  the  first  law  of  thermo¬ 
dynamics,  is  written  as  follows: 


d04  _  [  Qcnd,3-4  +  0cnd,4-5  +  Qohm,4  +  0mass,4] 

~dx  ~  W4 


(26) 


The  dimensionless  heat  conduction  rate  from  the  cathode  layer 
to  the  electrolyte  is  written  as: 


0cnd,4-5  — 


2  &4&5^cnd,  4— 5(^5  ^4) 

§5/<4  +  §4/<5 


where  the  effective  conduction  heat  transfer  area  is  Acnd  4_5  =  (1  - 

05)^Z. 


where  the  dimensionless  convective  heat  transfer  rate  from  the 
oxidant  channels  to  the  cathode  layer  is: 


0cnv,5-6  =  fr6ttch?ch£z(l  -  05 X #6  “  O5) 


(30) 


And  the  energy  change  due  to  mass  transfer  across  the  cathode 
layer  boundaries  is: 

~  ho2(Ho2i6  -  Hq2)5)  +  ho2(XN2/Xo2)(HN256  -  Hn2;5) 

Omass,5  —  - ^ ~ ^ -  (31) 

^ref^p,f  J  ref 

The  energy  balance  in  the  oxidant  channels  in  CV6  yields: 

d0  6  [— Qcnv,l-6  —  0cnv,5-6  +  Qmass,6] 


dr 


N6 


(32) 


where  the  energy  change  due  to  mass  transfer  across  the  oxidant 
channels  boundaries  is: 


Qmass,6  — 


no2,in(Ho2  ,air  f^02,6)  +  ^N2,in(^N2,air  Hn2,6) 


(27)  JV6  = 


^refCp,flref 

The  dimensionless  heat  capacity  of  CV6  is  given  by: 

Pair^ch^/ch^/6f'z(Xo2  Cv?o2  TXn2Cv,N2  ) 


PTg/TIrefCpf 


(33) 


(34) 


The  pumping  power  needed  to  supply  fuel  and  oxidant  to  the 
single  cell  are  calculated  based  on  the  assumption  that  the  pumping 
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process  is  adiabatic: 


WP  =  m2 


+  m6 


k2  -  1 
*6 


(k2-\)/k2 


«2T2[(sii^a)  -i 

p^  +  ApsV"6-1^6  ' 


kq-  1 


^6^6 


Pair 


V 


-  i 


(35) 


where  the  pressure  drop  in  the  fuel  channels  and  oxidant  channels 
are: 


2 f2Lyp2u22 

Ap2  =  — n - ’ 

uh,2 


Ap6 


2f6Lzp6u2 

Dh,6 


(36) 


The  average  mass  flow  rates  in  the  fuel  and  oxidant  channels  are 
given  by: 

m2  =  2  [^H2(”H2,in  +  r*H2,out)  +  ^h2o^h2o]  (37a) 

rn6  =  1  [Mo2(no2,in  +  no2,out)  +  MN2(%2,in+nN2,out)]  (37b) 


In  Eq.  (35),  Kj  is  the  isentropic  index  and  given  as  kj  =  cpj7cvj-, 
j  =  2,6.  Pj  is  the  specific  gas  constant  and  calculated  by  Rj  =  cpj  -  cVJ-. 
Note  that  the  subscript  j  here  indicates  the  gas  mixture  in  the  fuel 
and  oxidant  channels. 


3.  Electrochemical  model 


To  complete  the  model  of  the  single  SOFC,  the  thermal  model 
needs  to  be  coupled  with  the  electrochemical  counterpart  before 
the  SOFC  performance  can  be  evaluated. 

The  actual  operating  potential  of  the  single  SOFC,  V,  is  derived 
from  the  open  circuit  potential,  V0cp,  by  taking  into  account  vari¬ 
ous  kinds  of  losses,  i.e.,  the  activation  losses  and  the  concentration 
losses  in  the  anode  layer  and  cathode  layer,  pactJ  and  pconcj,  as  well 
as  the  ohmic  losses  in  all  the  solid  compartments  of  the  fuel  cell, 

Pohmj- 


V  —  Vocp  —  I  'y  ^PohmJ  +  P cone,  3  +  Pc  one,  5  +  Pact,  3  +  Pact,  5  I  » 


j=  1,3, 4, 5 


(38) 


where  the  open  circuit  potential,  V0cp*  is  estimated  by  the  differ¬ 
ence  of  Gibbs  free  energy  between  the  reactants  and  products  of 
the  electrochemical  reactions  described  by  Eqs.  (1)  and  (2).  Its  cal¬ 
culation  is  reported  by  Kim  et  al.  [18]  based  on  the  assumption  that 
both  reactions  occur  at  the  electrolyte  temperature,  T4: 


V0CP  = 


-{(Hh2o,4  -  Hh2,4  -  i/2Ho2,4)  -  r4(s£20,4  -  s°2>4  -  1/2S°2  4 

2 F 


The  total  potential  losses  in  each  electrode,  ptotaij»  is  defined  as: 

Ptotalj  —  Pact,j  +  PohmJ  +  Pconc,j>  J  —  3,  5  (40) 

The  composite  electrodes  of  the  single  SOFC  are  a  mixture  of 
electronic  and  ionic  conductor  (electrolyte  material).  In  the  present 
study,  the  electrolyte  material  in  the  composite  electrodes  is  mod¬ 
eled  as  regularly  spaced  corrugations,  which  protrude  from  the 
dense  electrolyte  surface,  as  shown  in  Fig.  2.  The  electronic  conduct¬ 
ing  material  is  spread  over  the  electrolyte  material.  The  remaining 
space  corresponds  to  parallel  tube-like  pores  for  the  transport  of 
gaseous  species  [21  ].  The  interface  of  the  ionic  conductor,  electronic 
conductor  and  reacting  gases  is  called  three-phase  boundary  (TPB), 
and  corresponds  to  the  sites  where  the  electrochemical  reaction 
could  take  place.  The  active  TPB  regions,  L'j  opt,  refer  to  the  TPB 
sites  where  the  electrochemical  reactions  actually  occur.  Accord¬ 
ing  to  Chan  and  Xia  [22],  the  electrochemical  reactions  take  place 


where  the  total  potential  losses  in  each  electrode,  ^otaij.  reach  a 
minimum. 

The  ohmic  loss,  pohmj-,  is  obtained  based  on  the  Ohm’s  law  as 
follows: 


Pohmj  —  Wj  >  J  —  1  >  4 


(41a) 


_  WW-WS'opt  '(ij-ij.opt) 
VohmJ  2Aj  optLyIz(l  -  0,- )<t4  +  LyLz(  1  -  1 .80,  )aj  ’ 


j  =  3,5  (41b) 


where  the  electrical  resistances  of  the  interconnect  and  electrolyte, 
/3i  and  fi4,  are  calculated  according  to  their  material  and  geometric 
features  as  follows  [1]: 


a 


bi _ h _ _ _ 

Iylzai  ( 1  +  nch  )LtLya-[  ( 1  +  nch  )LtLza \ 


(42a) 


P*  = 


u 

LyL^CT/\ 


(42b) 


where  exj  is  the  electrical  conductivity  of  each  control  volume. 

The  first  term  at  the  right  side  of  Eq.  (41b)  accounts  for  the 
ohmic  losses  in  the  active  TPB  regions  where  the  ionic  conduct¬ 
ing  resistance  is  dominant,  the  second  term  refers  to  the  remaining 
regions  of  the  electrodes  where  only  electronic  conducting  resis¬ 
tance  exists.  opt  refers  to  the  area  of  the  active  TPB  regions  in  each 
electrode,  and  it  is  calculated  with  A'j ?opt  =  (1  +  2U  opt0J/Dpj)LyLz. 
Here,  Dpj  stands  for  the  average  diameter  of  the  parallel  pores  in 
the  electrodes  and  it  is  obtained  with  Dpj  =DWJ0j/(l  -</>/).  DWJ  is 
the  width  of  a  corrugation  with  electrolyte  material  as  shown  in 
Fig.  2  and  its  value  is  given  by  Tanner  et  al.  [21  ]. 

When  electrochemical  reactions  take  place  rapidly,  the  con¬ 
centrations  of  reactants  at  the  reacting  region  drop,  which  result 
in  the  so-called  concentration  losses.  The  working  current  densi¬ 
ties  at  which  the  fuel  concentration  and  the  oxidant  concentration 
decrease  to  0  are  defined  as  the  anode  and  cathode  limiting  current 
densities,  ias  and  iCSt  respectively.  Based  on  the  TPB  pressure  corre¬ 
lation  given  by  Eq.  (21),  the  limiting  current  densities  are  written 
as: 


2FXH2pajrD303 
RT3  L3 

4FD505Pair  f _ Pair _ A 

y  Pair  —  ^02 Pair  J 


(43a) 

(43b) 


)  +  RT4  ln[PH2o/(PH2(po2)1/2)]) 


(39) 


The  concentration  overpotentials  in  the  anode  layer  and  cathode 
layer  are  given  by  Aguiar  et  al.  [19]  as: 


Pconc,c 


(44a) 


Pconc,a 


RT3  ,  /  Ph2Ph2o,tpb  \ 
2P  \Ph2,tpbPh2o  J 


(44b) 


The  activation  losses,  pactJ,  represent  the  energy  needed  to 
activate  electrochemical  reactions.  In  the  present  study  they  are 
estimated  by  the  Butler-Volmer  equation  [23]  as  follows: 


h  =  *0  J 


exp 


j  =  3,5 


-  exp 


(  (1  -  oic)nF 

V  RTJ 


(45) 
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In  Eq.  (45)  n  takes  the  value  of  2,  ac  stands  for  the  charge  transfer 
coefficient,  and  when  it  takes  the  typical  value  of  0.5,  the  explicit 
form  of  Eq.  (45)  for  the  activation  overpotentials  is  given  by: 

^=SSinh"1(2fc)’  J'  =  3’5  (46) 

where  z0j  stands  for  the  exchange  current  densities  and  are 
reported  by  Aguiar  et  al.  [19]  as  follows: 

kj = exp  (-]§■)  -  j=3’5  w 

where  Sj  and  Ej  are  the  pre-exponential  factor  and  the  activation 
energy,  respectively.  Their  values  are  reported  by  Aguiar  et  al.  [19]. 

In  Eq.  (46),  the  working  current  densities  at  both  electrodes,  z), 
are  defined  as  ij  =  I/A'j  opt. 

With  the  knowledge  of  actual  operating  potential  calculated  by 
Eq.  (38),  the  electrical  power  of  the  single  SOFC  is  obtained  with: 

WE  =  VI  (48) 

And  the  net  power  (available  for  utilization)  of  the  single  SOFC 
is  derived  from  its  electrical  power  by  subtracting  the  required 
pumping  power: 


WN  =  WE-  Wp 


(49) 


The  working  current  and  calculated  potentials  and  powers  from 
the  model  are  then  nondimensionalized  as  follows: 

/  =  jE  (50) 

*ref 

where  Jref  is  a  specified  reference  current  and  I  is  the  dimensionless 
current. 


Vj  = 


ref 


n 'l 

ref 


(51) 


where  Vref  =  (zhrefcp  fTref)//ref  stands  for  the  reference  voltage.  V 
and  fj  are  the  dimensionless  potential  and  potential  losses,  respec¬ 
tively.  Subscript  j  denotes  all  the  potentials  and  potential  losses 
present  in  the  fuel  cell. 

Wi 

wj  =  Tjf~’  J  =  E,P,N  (52) 

Wref 

where  Wref  is  the  reference  power  defined  as  Wref  =  rhrefCp  fTref.  W 
is  the  dimensionless  power. 


4.  Results  and  discussion 

In  an  effort  to  validate  the  model,  we  compared  the  electro¬ 
chemical  performance  simulated  with  the  current  model  with 
results  published  by  Aguiar  et  al.  [19].  Fig.  3  illustrates  such 
comparison.  It  can  be  seen  that  under  same  SOFC  geometric 
dimensions  and  operating  condition  (1073 1<),  the  general  shape  of 
two  polarization  curves  are  similar,  as  well  as  the  relative  magni¬ 
tude  between  each  kind  of  potential  losses.  The  biggest  difference 
observed  is  in  the  ohmic  losses.  This  can  be  explained  from  the  dif¬ 
ference  in  the  modeling  of  ohmic  losses.  The  main  difference  is  that 
in  the  present  model:  (i)  the  composite  electrodes  are  a  mixture 
of  electronic  and  ionic  conductor  with  different  conductivities; 
(ii)  the  charge-transfer  reaction  happens  at  the  active  TPB  region 
expanding  a  small  distance  from  the  electrode/electrolyte  inter¬ 
face;  (iii)  the  ionic  conducting  resistance  in  the  active  TPB  region  is 
accounted  for;  (iv)  ohmic  losses  of  the  interconnect  are  accounted 
for  in  the  current  model.  As  a  result,  the  ohmic  losses  in  the  present 
model  are  slightly  larger  than  the  one  in  Aguiar  et  al.  [  1 9].  However, 


Table  1 

Properties  and  constant  values  used  in  the  simulations. 


B  =  5.6697  x  1 0“8  W  irr 2  I<“4 

Vr  =  3.07  x  10-5  m3 

cP)i=0.5xl03Jkg-1I<-1 

N 

II 

N 

II 

O 

cn 

3 

cPj3  =0.595  x  103  J  kg-1  K"1 

a  =  0.5 

cp’4  =  0.606  x  103Jkg“1  K"1 

83  =6.54  x  1011  £2_1  irr2 

cP)5  =  0.573  xl03J  kg-1  K-1 

<55  =  2.35x1011^-1  m-2 

cP)f  =  1 4983  J  kg-1  K_1 

£i=0.1 

cv,h2  =  21.718J  mol-1  IC1 

£3  =  £5  =0.8 

cv,h2o  =  33.156J  moP1 1<-1 

Vo2  -  16.3  cm3 

cv,o2  =  26.56 J  mol”1  IC1 

Vn2  =  18.5  cm3 

cv,n2  =  24.36J  mol-1  K"1 

Vh2  -  6.12  cm3 

Dw,3  =  14  x  10“6  m 

Vh2o  =  13.1  cm3 

Dw, 5  =  14  x  10“6  m 

£ch  =  6.4xl0-2 

£3  =  140  x  103  J  mol-1 

£t  =  2.56  x  10_1 

£5  =  137  x  103J  moP1 

pi  =  8000  kg  irr3 

f=  96500  C  moP1 

p3  =  6870  kg  m'3 

/ref-  116.58  A 

p4  =  5900  kg  m3 

kf=k2  =  0.43  Wm-1  K"1 

p5  =  6570  kg  m~3 

kair  =  /<6  =  0.0717Wm-1  K1 

<j\  =  1.5  x  106  ^_1  m_1 

k\  =25  Wm-1  K-1 

cr3  =  8.0  x  104  ^_1  m_1 

k3  =  /<5  =  2  Wm1  K1 

cr4  =  33.4  x  103  exp(-10.3  x  10 3/r)^“1  m”1 

/<4  =  2  Witt1  K-1 

(j 5  =  8.4  x  103  ^_1  m_1 

nch  =  10 

03=0.4 

Pf  =  Pair  =  0.1MPa 

05=0.5 

pref  =  0.1  MPa 

03=9.5 

ft  =  8.314J  mol-1  K"1 

05  =  7.2 

lair.in  =  923.15  I< 

/Xf  =  2.07  x  10“5  Pas 

rfiin  =  923.15  K 

/xair  =  4.18  x  10“5  Pas 

rref  =  298.15 1< 

the  qualitative  behavior  of  both  electrochemical  models  is  very 
similar. 

The  values  of  the  physical  properties  utilized  to  obtain  the  sim¬ 
ulation  results  shown  in  this  section  with  the  present  model  are 
listed  in  Table  1.  The  inlet  temperature  and  pressure  of  both  fuel 
and  oxidant  are  assumed  as  923.1 5 1<  (Tf,Tair)  and  1  atm  (Pf,pair), 
respectively.  The  geometric  features  of  the  single  SOFC  are  shown 
in  Fig.  1,  which  include  the  external  dimensions  (L*,  Ly,  and  Lz), 
the  individual  layers  thicknesses  (Li),  the  width  (Lch)  and  shoulder 
thicknesses  (Lt)  of  fuel  and  oxidant  channels,  and  the  number  of 
channels  embedded  in  the  bipolar  plate  (rzch). 

As  discussed  in  the  previous  section,  the  thermal  model  of 
the  single  SOFC  consists  of  six  dimensionless  ordinary  differential 
equations  to  obtain  the  temperatures  of  each  CV  and  three  algebraic 
equations  to  obtain  the  partial  pressures  at  the  active  three-phase 
boundaries  (TPB)  regions,  i.e.,  Eqs.  (21a)-(21c).  The  present  study 
focuses  on  the  optimization  of  SOFC  internal  structure  and  oper¬ 
ating  conditions  at  steady  state  to  achieve  the  maximum  electrical 
power  output.  As  a  result,  only  the  steady  state  solutions  of  the 
model  are  considered,  which  are  obtained  by  making  the  time 
derivatives  of  the  temperature  for  each  control  volume  equal  to 
zero  and  then  solving  them  with  a  quasi-Newton  method  [24]. 
The  contribution  of  Qactconcj  was  assumed  negligible  in  Eqs.  (23) 
and  (29).  Activation  and  concentration  losses  still  affect  the  nom¬ 
inal  temperature  through  the  actual  operating  potential  Eq.  (38). 
The  obtained  solutions  are  then  used  in  the  electrochemical  model, 
with  which  the  electrochemical  performance  of  the  single  SOFC  is 
evaluated. 

Fig.  4  shows  the  way  to  determine  the  active  TPB  region  at 
the  anode,  §'3?opt,  which  in  this  model  corresponds  to  the  reaction 
layer  thickness,  §'3,  that  minimizes  the  total  potential  losses  at  the 
anode.  Three  total  anode  thicknesses  are  considered,  £3  =6.4x1 0-4 
(dotted  line),  £3  =  16  x  10-4  (dashed  line)  and  §3  =  160  x  10-4  (solid 
line)  (the  dotted  and  dashed  lines  overlap).  The  total  potential 
losses  at  the  anode  reach  a  minimum  at  §'3)0pt  =  3.52  x  lO-4,  3.2  x 
10-4  and  2.98  x  10-4  for  working  currents  of  I  =  0.172,  0.344 
and  0.516,  respectively,  regardless  of  the  total  anode  thickness 
§3.  The  optimal,  §'3)0pt,  is  insensitive  to  the  changes  in  the  anode 
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Fig.  3.  Comparison  of  polarization  curves  and  power  density  simulated  in  current  model  (left)  and  on  Aguiar  et  al.  [19]  (right). 
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Fig.  4.  The  total  potential  losses  at  the  anode  and  electrical  power  and  net  power 
as  functions  of  anode  reaction  region  thickness  under  various  working  currents  and 
anode  thickness. 


Fig.  5.  The  total  potential  losses  at  the  cathode  and  electrical  power  as  functions 
of  cathode  reaction  region  thickness  under  various  working  currents  and  cathode 
thicknesses. 


thickness.  Fig.  4  also  shows  that  the  value  of  where  the  total 
potential  losses  in  the  anode  reach  a  minimum  is  also  the  active 
anode  thickness  where  electrical  and  net  power  reach  a  maximum 
(e.g.,  ^3?0pt  =  2.98  x  10-4  for  a  working  current  of  I  =  0.516).  The 
electrical  power  and  net  power  curves  are  parallel,  since  for  the 
conditions  of  Fig.  4,  the  pumping  power  is  constant  for  each  current 
level. 

Fig.  5  shows  the  way  to  determine  the  active  TPB  region  at  the 
cathode,  §'5)0pt.  Similar  conclusions  to  the  anode  side  are  obtained: 
the  optimal  £'5,  is  practically  insensitive  to  the  changes  in  the  cath¬ 
ode  thickness;  the  optimal  % 5  where  the  total  potential  losses  in 
the  cathode  reach  a  minimum  is  also  the  optimal  where  electri¬ 
cal  and  net  power  reach  a  maximum.  When  the  working  currents 
I  are  0.172,  0.344  and  0.516,  the  active  TPB  regions  at  the  cathode, 
^5, opt’  are  5-19  x  10-4’  4.61  x  10-4  and  4.13  x  10-4,  respectively. 
Since  the  total  potential  losses  at  the  anode  and  cathode  are  com¬ 
puted  separately  and  do  not  influence  each  other,  %'3,opt  and  t 
are  independent  of  each  other. 


Fig.  6  shows  the  typical  temperature  profile  of  each  compart¬ 
ment  in  the  single  SOFC.  It  is  observed  that  as  the  current  increases, 
the  temperature  of  all  the  control  volumes  increase  as  well.  The 
reason  is  that  the  higher  the  current  is,  the  more  fuel  is  consumed 
in  the  electrochemical  reaction,  as  a  result,  more  heat  is  gener¬ 
ated  in  the  single  SOFC  which  causes  the  temperature  of  all  CVs 
to  increase.  Fig.  6  also  shows  that  there  is  significant  temperature 
spatial  gradient  across  each  compartment,  especially  the  temper¬ 
ature  difference  between  the  oxidant  channels  (6)  and  the  PEN  (3, 
4,  5)  (positive  electrode-electrolyte-negative  electrode  assembly) 
in  the  present  model.  This  is  due  to  the  poor  thermal  conductivity 
of  oxygen  compared  with  hydrogen  and  the  solid  materials  of  the 
SOFC.  The  temperature  spatial  gradient  has  a  direct  impact  on  the 
fuel  cell’s  performance,  therefore  the  performance  evaluation  (fea¬ 
tured  by  electrical  power  and  polarization  curves)  should  take  all 
the  internal  spatial  temperature  gradients  into  consideration. 

Fig.  7  shows  the  typical  polarization  curves  of  the  single  SOFC. 
As  the  current  increases,  the  temperatures  of  all  the  SOFC  com- 
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Fig.  6.  The  temperature  of  each  compartment  in  the  single  SOFC  as  functions  of 
current. 


Fig.  8.  The  working  current,  electrical  power  and  electrolyte  temperature  at  the 
nominal  operating  point  as  functions  of  air  stoichiometric  ratio. 


Fig.  7.  The  polarization  curves  of  the  single  SOFC. 
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Fig.  9.  The  electrolyte  temperature,  electrical  power  and  net  power  at  the  nominal 
operating  point  as  functions  of  fuel  utilization  factor. 


partments  increase.  As  a  result,  the  Gibbs  free  energy  difference 
between  the  reactants  and  products  of  the  electrochemical  reac¬ 
tions  decreases,  which  leads  to  the  drop  in  open  circuit  potential 
according  to  Eq.  (39).  In  addition,  all  types  of  overpotentials, 
including  the  activation  overpotential,  concentration  overpoten¬ 
tial  and  ohmic  overpotential,  increase  monotonically  as  the  current 
increases,  which,  combined  with  the  decrease  of  open  circuit  poten¬ 
tial,  results  in  the  drop  of  the  actual  working  potential,  according 
to  Eq.  (38).  Besides,  it  is  observed  in  Fig.  7  that  for  small  currents, 
the  activation  overpotential  is  dominant  among  all  the  potential 
losses.  With  the  increase  of  current,  the  ohmic  overpotential  grows 
rapidly. 

Fuel  cells  are  normally  designed  to  operate  at  a  cell  voltage 
between  0.6  and  0.7  V  [19]  and  specifically,  SOFC  systems  are 
presently  rated  at  a  cell  voltage  of  0.7  V  per  cell  [25].  In  this  study, 
the  nominal  operating  point  is  defined  as  the  operating  point  at 
which  the  actual  working  potential  is  0.7  V.  From  the  polarization 
curves  of  Fig.  7,  the  dimensionless  working  current  at  the  nominal 
operating  point  is  1.1,  and  from  Fig.  6,  the  dimensionless  tempera¬ 
ture  of  the  PEN  at  the  nominal  operating  point  is  about  3.72. 


The  effect  of  air  stoichiometric  ratio,  A,air,  on  the  performance  of 
the  single  SOFC  is  illustrated  in  Fig.  8.  It  is  observed  that  the  nominal 
current  and  electrical  power  decrease  as  the  air  stoichiometric  ratio 
increases,  which  illustrates  that  large  air  stoichiometric  ratios  do 
not  benefit  SOFC  performance.  Such  phenomenon  can  be  explained 
physically:  as  the  air  stoichiometric  ratio  increases,  more  air  with 
comparably  low  temperature  flows  through  the  oxidant  channels, 
which  leads  to  the  drop  in  SOFC  temperature,  as  shown  in  Fig.  8. 
As  the  temperature  of  the  single  SOFC  decreases,  the  total  potential 
losses  increase  rapidly  and  the  actual  working  potential  drops,  end¬ 
ing  up  with  small  current  and  low  electrical  power  at  the  nominal 
operating  point. 

Fig.  9  illustrates  the  significant  impact  that  the  fuel  utilization 
factor  has  on  the  performance  of  the  single  SOFC.  The  maxi¬ 
mum  dimensionless  nominal  electrical  power  reaches  0.765  when 
yf  =  0.88,  which  is  approximately  110%  higher  than  the  nominal 
electrical  power  of  0.362  when  yf  =  0.1.  Fig.  9  also  shows  that  very 
large  or  very  small  fuel  utilization  factor  does  not  benefit  the  SOFC 
performance.  This  effect  can  be  explained  physically  by  analyzing 
two  extremes:  (i)  for  large  fuel  utilization  factor,  the  potential  losses 
decrease  due  to  the  relatively  high  temperature  in  SOFC  (less  fuel 
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Fig.  10.  The  electrical  power  at  the  nominal  operating  point  as  a  function  of  anode 
reaction  region  thickness  for  various  electrolyte  thicknesses. 

with  comparably  low  temperature  at  inlet  flows  through  the  fuel 
cell),  as  shown  in  Fig.  9,  however,  due  to  the  low  concentration 
of  hydrogen  in  the  fuel  channels,  the  open  circuit  potential  also 
decreases  according  to  Eq.  (39),  which  leads  to  the  drop  of  actual 
working  potential  and  ends  up  with  small  nominal  current  and  elec¬ 
trical  power;  (ii)  for  small  fuel  utilization  factor,  the  open  circuit 
potential  increases  (due  to  the  high  concentration  of  hydrogen  in 
the  fuel  channels),  but  the  potential  losses  become  large  due  to  the 
low  temperature  in  the  fuel  cell  (more  fuel  with  comparably  low 
temperature  at  inlet  flows  through  the  single  SOFC),  so  the  actual 
working  potential  decreases  at  this  extreme  too,  which  leads  to  the 
low  nominal  current  and  electrical  power. 

Under  a  constrained  operating  condition  (fixed  fuel  utilization 
factor  and  air  stoichiometric  ratio),  the  optimization  problem  of 
the  PEN  internal  structure  consists  of  pursuing  maximum  nom¬ 
inal  current  and  electrical  power  by  optimizing  the  thickness  of 
the  anode,  electrolyte  and  cathode  layers  subject  to  a  fixed  total 
thickness  of  PEN,  i.e.,  L3  +  L4  + 15  =  constant.  Note  that,  for  most 
single  SOFCs  produced  currently,  the  minimum  thickness  of  the 
anode,  electrolyte  and  cathode  layers  are  50  p,m,  5  pum  and  50  pan, 
respectively,  due  to  the  limitation  of  manufacturing  techniques. 
The  analysis  and  discussion  about  the  PEN  internal  structure  in  this 
paper  is  based  on  a  more  extended  range,  and  the  conclusion  from 
it  could  be  used  as  reference  in  the  future.  The  thickness  ranges 
of  the  anode  layer,  electrolyte  and  cathode  layer  are  5-560  pan, 
5-30  am  and  5-560  am,  respectively,  and  L3  +L4  +  L5  =  570  am 
(?3+?4  +  ?5  =  182.4  xlO-4). 

Fig.  10  shows  that  under  a  fixed  operating  condition,  §'3)0pt,  i.e., 
the  optimal  £'3  in  which  the  nominal  electrical  power  reaches  a 
maximum,  is  robust  with  respect  to  changes  in  electrolyte  thick¬ 
ness:  £'3>opt  is  2.399  x  10-4, 2.415  x  10“4  and  2.431  x  10~4  when§4 
is  1.6  x  10-4, 4.8  x  10-4  and  8  x  10-4,  respectively.  With  the  same 
method  the  effect  of  electrolyte  thickness  on  £'5}0pt  is  obtained: 
£'5t0pt  is  2.822  x  10-4,  3.129  x  10“4  and  3.332  x  10"4  when  §4  is 
1.6  x  10-4, 4.8  x  10-4  and  8  x  10"4. 

Fig.  1 1  shows  how  and  $'5  opt  change  as  the  electrolyte 
thickness  varies.  It  is  shown  that  with  the  selected  thickness  range 
of  electrolyte,  the  active  TPB  region  at  the  cathode,  f'5j0pt.  is  larger 
than  the  counterpart  at  the  anode;  as  the  electrolyte  increases,  both 
^  3, opt  an(l  §5  opt  increase>  though  §'3)0pt  is  practicably  insensitive 
(i.e.,  robust)  with  respect  to  electrolyte  thickness  variation.  The 


Fig.  11.  The  active  reaction  region  thickness  at  the  anode  and  cathode  and  the  two- 
way  maximized  electrical  power  at  the  nominal  operating  point  as  functions  of 
electrolyte  thickness. 

results  of  Fig.  1 1  are  used  in  the  next  discussion  to  set  the  values  of 
£  3, opt  and  §5  opt-  The  two-way  maximized  electrical  power  at  the 
nominal  operating  point  is  also  reported  in  Fig.  11. 

Fig.  12  illustrates  the  impact  of  thickness  distribution  between 
two  electrodes  on  the  performance  of  the  single  SOFC,  when  the 
electrolyte  thickness  is  fixed.  For  different  electrolyte  thicknesses, 
the  nominal  electrical  power  increases  as  the  anode  increases  at 
first,  and  then  starts  to  drop  after  the  anode  thickness  increases 
to  an  optimal  value.  The  optimal  thickness  distribution  between 
two  electrodes  is  determined  by  the  trade-off  between  the  poten¬ 
tial  losses  at  the  anode  and  cathode:  (i)  for  small  anode  and  large 
cathode,  the  potential  losses  at  the  anode,  primarily  the  concentra¬ 
tion  losses,  would  decrease,  but  the  potential  losses  at  the  cathode 
would  increase,  mainly  the  ohmic  losses;  therefore,  the  total  poten¬ 
tial  losses  of  the  SOFC  would  increase,  which  leads  to  the  drop  of 
nominal  electrical  power;  (ii)  for  the  large  anode  and  small  cath¬ 
ode,  the  potential  losses  at  the  cathode  would  decrease,  especially 
the  ohmic  losses,  but  the  potential  losses  at  the  anode,  mainly  the 
concentration  losses,  would  increase.  In  this  way,  the  total  poten¬ 
tial  losses  of  the  SOFC  would  increase  at  this  extreme  too,  leading 
to  a  drop  of  nominal  electrical  power. 

Fig.  13  shows  the  optimal  thickness  distribution  between  the 
anode  and  cathode  as  a  function  of  electrolyte  thickness  under  a 
fixed  operating  condition.  With  the  selected  thickness  range  of  elec¬ 
trolyte,  §3)0pt  is  much  larger  than  f5t0pt,  which  suggests  that  under 
the  fixed  operating  condition,  an  anode-supported  SOFC  would  per¬ 
form  better  than  a  cathode-supported  SOFC.  Fig.  1 3  also  shows  that 
optimal  anode  and  cathode  thicknesses  are  essentially  unchanged 
by  variations  in  the  electrolyte  thickness.  As  the  electrolyte  thick¬ 
ness  increases,  the  optimal  anode  thickness  increases  slightly  and 
becomes  stable  when  the  electrolyte  is  larger  than  8  x  10-4.  The 
opposite  trend  is  observed  for  the  optimal  cathode  thickness:  as  the 
electrolyte  grows,  the  optimal  cathode  thickness  decreases  to  some 
extent  at  first  and  becomes  stable  when  the  electrolyte  is  larger 
than  8  x  10-4.  These  observations  corroborate  the  robustness  of 
the  optimal  internal  structure. 

Fig.  13  also  shows  the  important  impact  the  electrolyte  has  on 
the  performance  of  the  single  SOFC.  Under  the  fixed  operating  con¬ 
dition,  for  the  optimal  electrode  thicknesses  distribution  and  the 
optimal  active  reaction  region  thickness,  the  three-way  maximized 
nominal  electrical  power  drops  sharply  as  the  electrolyte  thick¬ 
ness  increases,  due  to  its  large  ionic-conducting  resistance  and  the 
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Fig.  12.  The  two-way  maximized  electrical  power  at  the  nominal  operating  point  as  a  function  of  anode  layer  thickness  under  various  electrolyte  thicknesses. 


Fig.  13.  The  optimal  thickness  distribution  between  the  anode  and  cathode  as  well 
as  the  three-way  maximized  electrical  power  at  the  nominal  operating  point  as 
functions  of  electrolyte  thickness. 


resulting  large  ohmic  losses.  Therefore,  the  improvement  of  the 
ionic-conductivity  of  the  electrolyte  and  the  decrease  of  its  thick¬ 
ness  are  key  issues  to  improve  SOFC  performance. 

So  far,  the  internal  structure  optimization  has  been  pursued  at 
fixed  operating  conditions  (i.e.,  yf  =  0.75,  Aair  =  7).  In  the  remain¬ 
ing  part  of  this  study,  the  assumption  of  constant  fuel  utilization 
factor  is  relaxed.  In  practice,  since  the  minimum  thickness  of  the 
electrolyte  is  determined  by  the  manufacturing  techniques,  the 
optimization  problem  is  simplified  as  follows:  for  a  fixed  electrolyte 
thickness  and  given  total  PEN  thickness,  the  optimal  thickness  dis¬ 
tribution  between  the  anode  and  cathode  is  pursued  to  achieve  the 
maximum  nominal  electrical  power.  The  four  degrees  of  freedom 
are  then  the  thickness  of  the  active  reaction  layer  region  in  the 
anode  and  cathode,  the  anode  and  cathode  thickness  distribution, 
and  the  fuel  utilization  factor. 

Fig.  14  shows  how  the  nominal  electrical  power  changes  with 
the  thickness  of  anode  reaction  region  under  various  fuel  utilization 
factors.  §3,0pt»  the  optimal  %3  where  the  nominal  electrical  power 
reaches  maximum,  changes  slightly  for  different  fuel  utilization  fac¬ 
tors:  £'3>0pt  is  3.11  x  10-4,  2.59  x  10~4  and  2.39  x  10~4  for  yf  =  0.1, 
0.5  and  0.9,  respectively.  Similarly,  for  the  cathode  side,  ^opt  is 
3.18  x  10-4,  3.22  x  10-4  and  3.26  x  10~4  for  yf  =  0.1,  0.5  and  0.9, 
respectively. 
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Fig.  14.  The  electrical  power  at  the  nominal  operating  point  as  a  function  of  anode 
reaction  region  thickness  under  various  fuel  utilization  factors. 


Fig.  16.  The  two-way  maximized  electrical  power  at  the  nominal  operating  point 
as  a  function  of  fuel  utilization  factor  for  various  anode  layer  thicknesses. 


Fig.  15  shows  how  tj'3  opt  and  opt  change  as  the  fuel  utiliza¬ 
tion  factor  varies:  as  the  fuel  utilization  factor  increases,  the  active 
reaction  region  at  the  anode  decreases,  whereas  the  active  reaction 
region  at  the  cathode  increases  slightly.  The  results  in  Fig.  15  are 
used  in  the  next  discussion  to  set  the  values  of  §'3j0pt  and  $'5  opt  as 
functions  of  the  fuel  utilization  factor. 

Fig.  1 6  shows  how  the  two-way  maximized  nominal  electrical 
power  (e.g.,  corresponding  to  £'3>opt  and  t)  changes  with  fuel 
utilization  factor  for  different  thickness  distributions  between  the 
two  electrodes.  The  optimal  fuel  utilization  factor,  y fopt,  are  0.844, 
0.863  and  0.878  when  §3  is  1.6  x  10"4,  64  x  10“4  and  160  x  10-4, 
respectively.  The  resulting  relationship  between  the  optimal  fuel 
utilization  factor  and  the  thickness  distribution  between  the  two 
electrodes  is  shown  in  Fig.  17.  It  is  observed  that  yfi0pt,  the  opti¬ 
mal  fuel  utilization  factor,  increases  as  the  anode  layer  thickness 
increases. 

Since  ij'3j0pt  and  §5  t  are  functions  of  the  fuel  utilization  fac¬ 
tor  (Fig.  15)  and  the  optimal  fuel  utilization  factor  is  a  function 
of  the  thickness  distribution  between  the  two  electrodes  (Fig.  17), 
then  £3  opt  and  opt  are  functions  of  the  thickness  distribution 


between  two  electrodes  under  optimal  fuel  utilization  factor.  This 
relationship  is  shown  in  Fig.  18. 

The  combination  of  Figs.  1 7  and  1 8  shows  that  as  the  anode  layer 
thickness  increases,  the  optimal  fuel  utilization  factor  increases  as 
well,  and  accordingly,  §'3>opt,  the  active  TPB  region  at  the  anode, 
decreases  slightly,  and,  $'5j0pt,  the  active  TPB  region  at  cathode 
increases  slightly. 

The  three-way  maximized  nominal  electrical  power  as  a  func¬ 
tion  of  the  thickness  distribution  between  the  two  electrodes  is 
shown  in  Fig.  17.  With  the  selected  thickness  range  of  anode,  the 
k^E,n,mmm  increases  as  the  anode  thickness  increases,  which  means 
that  under  flexible  operating  conditions,  the  anode-supported  SOFC 
exhibits  better  performance  than  the  cathode-supported  SOFC, 
when  all  other  parameters  are  the  same. 

It  is  important  to  note  that  in  the  present  study,  the  properties 
were  evaluated  at  a  representative  temperature  and  treated  as  con¬ 
stants  (see  Table  1 ).  The  properties  dependence  on  temperature  can 
be  accounted  by  introducing  appropriate  empirical  correlations.  To 
quantify  the  sensitivity  with  respect  to  temperature  of  the  param¬ 
eters  shown  in  Table  1  and  their  effect  on  the  simulation  results, 


Fig.  15.  The  active  reaction  region  thicknesses  at  the  anode  and  cathode  as  functions  Fig-  17.  The  optimal  fuel  utilization  factor  and  three-way  maximized  electrical 

of  fuel  utilization  factor.  power  at  the  nominal  operating  point  as  functions  of  anode  layer  thickness. 
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Fig.  18.  The  active  reaction  region  thicknesses  at  the  anode  and  cathode  as  functions 
of  anode  layer  thickness. 

an  operating  range  of  1000-1143  K  was  selected.  In  that  range,  the 
heat  capacities  vary  by  4%,  the  thermal  conductivity  of  fuel  and  air 
vary  by  10.8%  and  their  viscosity  by  9%.  The  corresponding  impact 
in  the  active  TPB  thickness  at  both  electrodes  is  less  than  1%.  The 
variation  of  each  control  volume  temperature  is  less  than  0.5%. 
The  activation  overpotential  displays  a  variation  of  4.8%,  the  ohmic 
overpotenial  varies  by  2.5%  and  the  concentration  overpotential  by 
2.2%. 

The  temperature  profile  of  the  SOFC,  the  polarization  curves 
and  the  active  TPB  at  both  electrodes  are  the  foundation  for  the 
optimization  of  internal  configuration  and  operating  condition.  The 
above  discussion  shows  that  the  maximum  variation  observed  was 
less  than  5%.  This  leads  to  conclude  that  assuming  constant  ther¬ 
mal  properties  in  the  given  operating  range  is  acceptable.  This  is 
consistent  with  a  comment  by  Bove  and  Ubertini  [11]. 

5.  Conclusions 

In  this  paper  several  important  parameters  related  to  the  SOFC 
performance  are  identified:  the  internal  structure  of  the  positive 
electrode-electrolyte-negative  electrode  assembly  (PEN),  which 
includes  the  active  TPB  regions  at  the  anode  and  cathode,  and  the 
operating  condition  (fuel  utilization  factor  and  air  stoichiometric 
ratio).  The  analysis  of  their  effects  on  the  SOFC  performance  is  con¬ 
ducted  based  on  a  numerical  model  which  accounts  for  both  the 
thermal  and  the  electrochemical  aspects  of  the  single  SOFC.  The 
model  takes  into  account  relevant  aspects  concerning  fluid  flow, 
mass  balance,  heat  and  mass  transfer,  energy  conservation,  elec¬ 
trochemical  reactions,  and  takes  spatial  temperature  and  pressure 
gradients  into  consideration. 

The  optimization  of  PEN  internal  structure  includes  the  opti¬ 
mal  thickness  allocation  between  the  anode,  the  electrolyte  and 
the  cathode  so  that  several  effects  such  as  the  open  circuit  poten¬ 
tial  and  all  kinds  of  losses  are  appropriately  balanced  to  achieve 
maximum  nominal  electrical  power.  The  procedure  also  addresses 
the  determination  of  optimal  active  TPB  regions  in  the  electrodes, 
which  is  based  on  the  trade-offs  between  all  the  potential  losses  in 
the  electrodes.  The  active  TPB  thickness  at  the  anode  obtained  with 
the  present  model  is  in  the  range  7.32-1 1  pan,  and  the  thickness  of 


the  active  TPB  at  the  cathode  is  in  the  range  8.82-16.2  pan.  Those 
values  are  consistent  with  the  reaction  zone  penetration  into  the 
electrode  reported  in  Janardhanan  and  Deutschmann  [12]  and  Nam 
and  Jeon  [26]. 

The  optimization  of  the  fuel  utilization  factor  and  air  stoichio¬ 
metric  ratio  is  also  based  on  the  balance  between  the  effects  of  SOFC 
operating  temperature  and  the  fuel  and  oxidant  concentration.  The 
maxima  found  are  sharp  and  therefore,  the  corresponding  optimal 
geometric  and  operating  parameters  could  be  used  as  preliminary 
rules  in  practical  SOFC  design. 

The  present  study  shows  that  trade-offs  not  only  exist,  but 
also  serve  as  the  basis  on  which  the  optimization  of  SOFC  inter¬ 
nal  structure  and  operating  condition  is  implemented.  In  practice, 
such  trade-offs  are  different  from  case  to  case,  therefore,  they 
must  be  pursued  based  on  models  that  correspond  to  real  appli¬ 
cations.  However,  the  methodology  of  modeling  and  optimization 
presented  in  this  paper  is  generally  applicable  and  can  be  extended 
to  the  analysis  of  large  and  more  complex  SOFC  systems. 
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